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We prove in full generality the thermal operator representation for Matsubara sums in a rela- 
tivistic field theory of scalar and fermionic particles. It states that the full result of performing the 
Matsubara sum associated to any given Feynman graph, in the imaginary-time formalism of finite- 
temperature field theory, can be directly obtained from its corresponding zero-temperature energy 
integral, by means of a simple linear operator, which is independent of the external Euclidean 
energies and whose form depends solely on the topology of the graph. 
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I. INTRODUCTION 



Relativistic quantum field theory is the mathematical formalism that allows us to describe the interactions among 
elementary particles, according to the principles of quantum mechanics and special relativity. When considering 
the physics of very dense or very hot plasmas, such as the early universe or the quark-gluon plasma that should 
' be produced in heavy ion collisions, it becomes necessary to formulate the relevant questions in terms of thermally 
averaged quantities, according to statistical mechanics. Most of the physical information, for systems both in and 
slightly out of thermal equilibrium, can be obtained from the study of the so-called thermal Green functions, defined 
as thermal expectation values of a time-ordered product of field operators. One fruitful and widely used approach to 
study these Green functions is perturbation theory, where these are computed by means of a systematic expansion in 
terms of Feynman diagrams |2|. 

The simplest diagrammatic analysis is obtained in the so-called imaginary-time formalism, in which the diagrams 
have the same topology and are computed according to basically the same Feynman rules as in the zero-temperature 
, Euclidean theory, except for one very important difference: the energy component of both the external 4-momenta and 
the internal 4-momenta carried by the propagators is quantized (in slightly different ways, according to the nature 
bosonic or fermionic- of the associated particle), that is, it must be a Matsubara frequency 1]. Accordingly, the 
calculation of a loop diagram in the imaginary-time formalism of quantum field theory at finite temperature necessarily 
involves sums over internal Matsubara frequencies 0, an operation that we shall generically call the Matsubara sum 
associated with the graph. Although this sum can be computed in a number of ways, usually in a systematic fashion, 



such computations can become quite tedious for higher loop diagrams 0, E(- Another related difficulty of the the 
imaginary-time formalism is the separation of the vacuum contribution from the finite temperature corrections, since 
the naive zero-temperature limit of any given loop graph is of the indeterminate type • oo [j| . 
k> \ It was discovered some time ago|^ that, for some particular classes of diagrams in a scalar theory, the full result 
of performing the Matsubara sum associated with a Feynman graph could be completely determined from its zero- 
temperature counterpart, by means of a simple linear operator that depends on the topology of the diagram but is 
independent of the Euclidean energies carried by its external lines. Due to the simple and well-defined structure of 
this operator, dubbed the thermal operator, it was conjectured that this result should hold for all Feynman graphs. 

In the recent paper p| it has been partially shown that this conjecture is actually correct. The authors of reference 
P| have rescued from oblivion an old, systematic and very elegant method to perform the Matsubara sum associated 
with an arbitrary graph, due to M. GaudinQ, and have used it, as an illustration of the power of Gaudin's method, 
to show that the main part of our conjecture follows naturally from it. The proof given in reference [j| does address 
the relationship between the structure of the full thermal result and the zero-temperature result, but stops short 
of explicitly constructing the thermal operator that relates them. Also, it leaves undiscussed conjectures 2 and 3 
presented in reference [(j. 

In this letter we restate all the conjectures of referenceQ as full-blown theorems concerning general properties of 
the Feynman graphs in the imaginary-time formalism of thermal field theory, extending their validity to theories 
containing Dirac fields as well. Following reference [j| , we make full use of Gaudin's method to prove these theorems, 
although in a slightly different and more explicit form. 
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The structure of the rest of the paper is as follows: In section [D] we present, in the form of theorems, the general 
properties of the finite-temperature imaginary-time Feynman graphs that were previously presented as conjectures, in 
terms of the thermal operator, suitably extended to incorporate possible fermionic lines. In section ITTll we summarize 
Gaudin's method to perform the sums over Matsubara frequencies, which will be central to the proofs of the theorems, 
presented in sections ITVl fscalars) and IVl ( fermions) . Our conclusions are presented in section IVTl 



II. THE THERMAL OPERATOR 

In a scalar field theory, the mathematical expression corresponding to an amputated graph with V vertices, / 
internal lines, and external 4-momenta P a = (p Q , p Q ) has the form 

^ / [n^rcr nW^) (kt ,)] D(p,E,n (i) 

a J . =1 (2tt) 2E t v=1 

where A represents the coupling constant and 5* is the symmetry factor of the graph; is the spatial 3-momentum 
of the i-th internal line and Ei := (k 2 + mf) 1 ^ 2 is its associated kinematic energy; k„ denotes the total 3-momentum 
entering vertex v; the unsubscripted symbols p and E denote, respectively, the full set of Euclidean external and 
kinematic internal energies, p :— {pi,p2, ■ ■ ■ ,Pv} (with ^2i =1 Pi — 0) and E :— {E\, E2, ■ ■ . , Ej}\ and T is the 
temperature. The delta functions ensure conservation of spatial 3-momentum at each vertex, so that the integration 
measure reduces essentially to an integration over the 3-momcnta of the L = I— V + 1 independent loops. In the finite 
temperature Euclidean formalism all scalar lines, external and internal, carry discrete Euclidean energies which are 
integer multiples of 2-7rT. Each internal line has an associated Matsubara frequency, denoted by fcj = Lu ni := 2-nTrii. 
The D-function is given by the normalized Matsubara sum 

1 

D(p, E, T) = 1E T L II &(h,Ei)8(p, k), (2) 

{rii} i=l 

where 



7b 



II 2 ^' ( 3 ) 



i=l 



L is the number of independent loops in the graph, and A(ki,Ei) is the scalar propagator associated with the i-th 
internal line, with 

w^wh?' (4) 

The sums over each rii (i — 1, . . . , I) run from —00 to +00. The (^-function, with k — {fci, . . . , ki}, is a generalized 
Kronecker delta which ensures conservation of energy at each vertex. The topology of the diagram is totally contained 
in this generalized delta. 

In a theory containing both fermions and scalars, the structure of Q is basically unchanged, except for extra 
spin indices carried by the external fermionic lines and a possible extra sign associated with fermionic loops. The D- 
function is still given by © , except that each fermionic line carries a Matsubara frequency consistent with antiperiodic 
boundary conditions on the fermionic fields, 

k = u) n := 27rT(n+ 1/2), (5) 

entering through the fermionic propagator, which now has a matrix structure and depends explicitly on the spatial 
momentum k, 

x /1 1 \ to + ikf - k ■ 7 

A(fc ' k) = — WTeI — • (6) 

where the 7^ = (7 ,7) are the usual Minkowski space gamma matrices, obeying the standard anticommutation 
relations, {7^,7^} = 2rf v . 
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In order to state our general theorems in a way that is applicable to a general theory containing both scalars and 
fermions, it will be convenient to define the following generalized signed thermal occupation number function, which 
takes into account the statistics corresponding to both kinds of internal lines: 

J n(E), for a scalar line; 

^ ' ' \ —h(E), for a fermionic line, ^ ' 

where n(E) — (e@ E — l) 1 and n(E) — (e' 5 ' 6 + l) 1 are, respectively, the Bose-Einstein and Fermi-Dirac thermal 
occupation factors, for the case of vanishing chemical potential. 

Additionally, the theorems will be formulated in terms of a simple reflection operator on the space of functions of 
several variables, 

S E f(E,x) := f(-E,x), (8) 

where E stands for any one variable and x for all the others. 
The main result presented in this paper is enunciated in the next theorem. It basically states that the "energy 
part" of any Feynman graph in the finite-temperature imaginary-time formalism, represented here in terms of the 
D-function introduced in QJ, can be obtained directly from the corresponding zero-temperature energy integral. 

Theorem 1. [Thermal Operator Representation] The D-function defined in (J2J for an amputated Feynman 
graph can be expressed in the form 

D(p,E,T) = 6(E,T)D Q (u,E) , (9) 

uj—p 

where Dq(u>,E) is the D-function of the Euclidean zero-temperature graph and (D(E,T), the thermal operator, is the 
following linear operator: 

i 

6(E, T) := 1+ + St) + N h N h (1 + 4)(1 + 4) 

L ^ l2) (10) 

+ ■•■+ E' LK(i+4)- 

Here Ni = N(Ei), where N(E) is the generalized signed thermal occupation factor defined in Q; Si :— iSe ; , where 
Se is the reflection operator defined in (JSJ); the indices ii,i2, ■ ■ ■ run from 1 to I (the number of internal propagators) 
and the symbol (ii,.. ■ ,ik) stands for an unordered k-tuple with no repeated indices, representing a particular set of 
internal lines. The primes on the summation symbols imply that certain tuples (ii,...,ik) are to be excluded from the 
sums: those such that if we snip all the lines ii, . . . , ik then the graph becomes disconnected. 

Note that the operator (D(E, T) contains products of at most L thermal occupation factors N(Ei), since for a L-loop 
graph the maximum number of lines that can be snipped without disconnecting the graph is precisely L. This generic 
feature of the thermal graph in the imaginary-time formalism is of course well known. However, based on a general 
property of the zero-temperature D-function, formulated in theorem [3] below, it is also possible to use a modified 
thermal operator, which has a simpler algebraic form: 

Theorem 2. [Simpler form of the Thermal Operator] When acting on the zero-temperature D-function, 
Dq(p,E), the thermal operator 0(E,T) can be replaced by the the simpler 

i 

G ic (E,T) = l[[l + N i (l + S i )]. (11) 

Note that the operator 0*(E,T) in 111|) can be expanded as in (|1(J|> . with the only difference that the summation 
symbols will carry no primes, that is, all tuples . . . ,ik) (1 < k < I) will be allowed in the sum. Clearly, forms 11U|) 
and (|llf) of the thermal operator will be equivalent if we can prove that tuples associated with disconnected graphs 
(the ones excluded from the summations in (|TU|> ) give rise to operators that produce a vanishing contribution to the 
D-function in 10. This is the content of our last theorem: 
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Theorem 3. [Cut sets do not contribute] The zero-temperature D-function, Dq(lu,E), is annihilated by the 
operators 

MO):= n( 1+ 4), (12) 

iiec 

where C stands for a cut set of the graph, that is, any set of indices i±, . . . ,i k such that the graph becomes disconnected 
if the corresponding lines are snipped. 

The concept of cut set, as used here, bears no connection to the concepts of cut and cut diagrams as they are usually 
understood in diagrammatic quantum field theory. Cut sets are determined solely by the topology of the diagram, 
and have no further mathematical or physical meaning. 



III. GAUDIN METHOD 



Here we present a summary of Gaudin's method and its use in the computation of the Matsubara D-function, in 
the purely scalar case. For all the details see references and [|. The changes in the presence of fermionic lines are 
commented upon at the end of this section. 

Gaudin's method to perform the Matsubara sum associated to a given graph is based on two main ideas. The 
first is to make use of the spectral representation of the propagator, which puts the dependence on the Matsubara 
frequency linearly in the denominator. For the scalar case it reads 

» / ^ f°° dk° p (k°, k) , 

with the spectral function given by 

p(*°,k) = 2 7 re(fc ( V(fc 2 - m 2 ) (14) 

= ^r[S(k°-E k )-5(k°+E k )}. (15) 

e(k°) is the sign of k°. 

Once all the propagators have been represented by means of i|13fl . the Matsubara sum to be performed takes the 
form 

ri £n^br^> fc) ' (16) 

{iJi}i=l 1 

where the fc, = uj ni are the Matsubara frequencies and the k® are external real variables. The generalized Kronecker 
delta 8{p, k) enforces V — 1 independent linear relations satisfied by the Matsubara frequencies k, also involving the 
external Euclidean energies p, which we shall write (following reference 0) as 

R v (p,k) = 0, for v = 1,...,V- 1. (17) 

This system of linear equations allows us to solve for V — 1 of the I Matsubara frequencies in terms of a set of 
L = I — V + 1 independent ones. In general, there will be several distinct ways of choosing this set of independent 
Matsubara frequencies. As shown by Gaudin, there is a one-to-one correspondence between the collection of all possible 
sets of independent Matsubara frequencies and the set of all trees associated to the given (connected) diagram T. 

A tree is a set of lines of T joining all vertices and making a connected graph with no loops. Every tree T will 
contain V — 1 lines and its complement T (the set of lines which do not belong to T) will have L lines. The Matsubara 
frequencies corresponding to the L lines in T, denoted by ki , will constitute a set of independent Matsubara frequencies 
in terms of which the system Ijl7|l can be solved. The Matsubara frequencies associated with the lines of the tree, kj , 
with j s T , will be linear combinations of the independent Matsubara frequencies and the external Euclidean energy 
variables, 

kj^njfakt), jeT,lef. (18) 



FIG. 1: The trees of a simple two-loop graph, shown as dark lines. 



As a simple example, in figure ^ we show the three possible trees for the two-loop two- vertex graph shown. In this 
case, each tree T is composed by a single internal line (heavy line), whose Matsubara frequency can be expressed, after 
using energy conservation at one of the vertices, in terms of the two independent Matsubara frequencies associated 
with the (thin) lines that do not belong to the tree (these are the lines in T) and the external energy p. For instance, 
for the first tree we have k$ = p — k\ — fc 2 , etc. 

Gaudin's main insight is the following identity for the rational function appearing in : 

1 I 1 ! 

EIlj.o jr. S (P> k ) = ^2 IT 1,0 _ IcytJZ. EII j.o_ jr ■ ( 19 ) 

In our example in figure ^ this identity takes the form 

\p S(p-ki-k 2 - fc 3 ) 1 

ni j^ n3 (A° - ikx) (fc° - ik 2 ) (fc 3 ° - ik s ) ~ (fc? + k\ + k\ - ip) 

X (fc?-^l)(fc 2 °-^2) + (fc? - ifci) (fc§ - ifc 3 ) + „5 3 (fc 2 °- *k 2 )(k»- l k 3 )} ■ 

So, using Gaudin's decomposition (|19|l . for any given graph the full Matsubara sum is split into a number of other 
sums, one for each tree T, where each one of these sums is actually simply a product of independent sums over each 
of the Matsubara frequencies in T. Each independent sum has to be regulated. Gaudin assigns to each internal line 
a regulator e ikiTi , where the time r, is taken to zero at the end. The sum that results is well defined, 

T E^-=^(^V ?T S (20) 
m 1 

where n(k°) is the Bose-Einstein occupation factor. TJ is some linear combination of the times Tj associated to the 
lines that belong to the loop defined by I £ T, and e/ is the sign of T/. Only e; matters when the regulators are 
removed. 

Consequently, Gaudin's result for the D-function in the purely scalar case is 



D(p, E,T) = f[ r dkt P {kl Ek) { II uQ^nr II ( ^ k ") I 



(21) 



where = &f(p, —ik®) and p is the reduced spectral function, 



p(k°, E k ) = 2 -^p (fc°, k) = S(k° ~ E k ) - S(k a + E k ). (22) 



Gaudin's method applies equally well in the case the graph contains fermionic lines. As for the scalar propagator, 
the fermionic propagator © admits a spectral representation, 

A^kl.ffifM, ( 23 ) 
J^oo 2ir k - tuj n 

with 



p (k , k) = 27T£(fc ) + m) <5(/c 2 - m 2 ), 



(24) 
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where the metric signature is that of Minkowski space, ft = k^j^ = /co7° — k • 7 and k 2 = k 2 — k 2 . Note that 
the spectral representation of the fermionic propagator |JBJ "hides" the dependence on the Matsubara frequency k 
appearing in its numerator, so that the only difference with the scalar case is the extra matrix structure (ft + m). 

Gaudin's method to perform the Matsubara sum in the form through the decomposition in terms of trees, 

makes no reference to the bosonic or fermionic nature of the frequencies to be summed over. If an independent 
frequency ki (I E T) is fermionic, we shall need the sum 



T Edr-r = -^^V" T, ; (25) 



fc° - iki 



which can be obtained from the bosonic analog by writing fc ; ° — iki = (k,® — z7rT) — iki. h(k°) is the Fermi-Dirac 
function, 

*(*°) = (26) 
We see that the fermionic result (|25[) can be simply obtained from its bosonic analog H2U[) through the replacement 

n — > —n. 

Therefore, the result for the D- function in the presence of fermionic propagators will be identical in form to i|21[l . 
except that for the fermionic lines the reduced spectral function p will have an extra matrix factor, and the thermal 
factor ein(eik\ ?) will be replaced by — e;ri(eifc ; ). 

IV. PROOFS OF THE THEOREMS FOR THE SCALAR CASE 

In order to prove the theorems enunciated in section [HJ we shall need to identify the vacuum (T = 0) limit of 
Gaudin's result l|21|) for the Matsubara D-function D(p, E,T). The temperature dependence is solely contained in 
the thermal factors n(ezfc ; °). Following reference 0, we use the identity 

n(k°) = -9(-k°) + e(k°)n(\k°\), (27) 

from which it follows that 

emW) - -eie{-eikf)+e(kf)n(\kf\). 
Writing the reduced spectral function defined in l|22ll in terms of the reflection operator Se as 

p(k°,E)= (l - S E ) S(k° - E), (28) 
and using the notation Si = <Se, , we can perform the integrations over each of the variables kf (I G T) in Ij21|l as 
drfpik^EOemie^fik?) = (l-Si) ^ dkfS(kf - £,)«■!»(*!*?)/(*?) 

^ ' J — OO 

'l - S;) em(eiEi)f(Ei) 
= (l-Si) [— ej0(— ej.Ej) + e(Ei)n(\Ei\)] f{E{), 

where Ei is for the moment considered as an arbitrary real variable. We now apply the reflection operator Si explicitly 
and then use the fact that Ei is actually a positive quantity. The vacuum part is given by 



l-Si)[-€ie{-€iEi)f{Ei)} = -ei6{-eiEi)f{E l )+e l B{e l E l )f{-Ei) 

= f(-eiEt), 
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whereas the thermal part is given by 

(l-Si) mMEtDfiEt)} = e(Ei)n(\E l \)f(Ei) - e(-Ei)n(\-Ei\)f(-Ei) 

= ni f(Ei)+nif(-Ei) 

= m (l + £,)/(£,) 



(l + Si) /(-e,£,), 



since ^1 + «S;^ f(Ei) = (l + S^j f{—E{), and where we have denoted ni = n(E{). 
Therefore, 



D(p,E,T) = J2H + (l + S^Dj&E) 



T ; 6 r 



(29) 



where D$ (p, 7?) is the contribution to the vacuum D-function associated with the tree T: 



/oo 
dp^ik^E 



1 



]> k°-inJ{pMEi) 



( T Ej-mjip^EiY 



(30) 



Equations Ij29(l and (|30[1 are our starting points for the proofs of our three theorems. To start with, we notice that 
the function Dj(p, E) is annihilated by each of the operators (1 + Sj) with j E T, due to the identity 



1 + S E ) ( 1 - S E ) = 1 - S% = 0. 



(31) 



This allows us to extend the index of the product Y[ m <|29JI to all possible values, for any tree T: 



1 

[] [l + ni (l + £i)]D%(p, E) = [] [l + n 3 (l + 4)] Dl{p, E) 



3 = 1 



In this form, the operator acting on Dq (p, E) becomes T- independent, which allows us to move the sum over all trees 
in through it: 



D(p,E,T) = n[l + ^(l+4)]E^o T (P,£) 

3 = 1 T 

= 6*(E,T)D (p,E), 

and this proves theorem 2. 

Theorem 3 is readily proven by noticing that if C is a cut set of the graph T, then the operator 



A(C) = J] (l + ^ 



kec 



will contain at least one factor 



(i+4v) 



with jV £ T for every tree T. As we have already pointed out, this factor 

will annihilate the corresponding D^(p, E), and therefore Do(p, E) will be annihilated by A(C). This proves theorem 
3. 

Finally, as explained in section [HI theorem 1 follows directly from theorems 2 and 3. 
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V. EXTENSION TO FERMIONS 



As it was shown in section Hill in the presence of fermionic lines the Matsubara D-function has essentially the same 
structure as in the scalar case. Now we use the identity 

n(fc°) = 0(-fc°) + e(fc°)n(|fc o |), (32) 

which has the same contents as the well-known h(—E) = 1 — h(E). We note that both identities (|27() and (|32|l can 
be written in terms of the generalized occupation number function defined by JJJ as 

N(k°) = -e(-k°) + e(k°)N(\k°\), (33) 

and therefore the manipulation of section HVI are valid in general, with rii replaced by iVj. 

The only other difference with the scalar case is the spectral function. The reduced spectral function for fcrmions 

is 

p(k°, k, E k ) = (k , k) = (# + to) [6(k° - E k ) - S(k° + E k )] = (l - S Ek ) S(k° - E k ) {#+ m) . (34) 

The propagator for a fermionic line i will be represented in the form 

1 „ f°° 1 

A(d)i,kj) = ——(l — Si) / dk J{k i -E i )$ + m)-3 r -r^, (35) 

J-oo K i ~ luJ i 

which has the necessary structure for the derivations presented in the previous section to hold. Hence, in the presence 
of fermionic lines the Matsubara D-function will still have the form (|29|l . apart from the change n — > —h to account 
for the Fermi-Dirac statistics for fermions and some extra structure in the numerator of the right-hand side of i|30[l . 
Therefore, the proofs presented at the end of the previous section are unaffected. 



VI. CONCLUSIONS 



In this paper we have proven rigorously, extending the approach presented in reference |5J] , a very general property of 
Feynman diagrams in the imaginary-time formalism for finite-temperature relativistic field theories (including scalar 
and fermionic fields), previously put forward as a conjecture in reference 6J. This property states that the full result 
of performing the Matsubara sum associated to any given Feynman graph can be obtained from its zero-temperature 
counterpart, by means of a simple linear operator, given by l|10|) . whose form depends solely on the topology of the 
graph. 

The thermal operator (|1(J|) has the important feature of being independent of the discrete Euclidean energies carried 
by the external lines of the graph. It follows from this that all issues related to analytic continuations of imaginary-time 
formalism Green's functions ja, El 0, EI can be completely settled at the zero-temperature level. The implications 
of this fact as well as the connections of the thermal operator representation with the real-time formalism are under 
investigation and will be presented elsewhere. 

We have not studied the validity of the thermal operator representation in the case of gauge theories. It would 
be interesting to determine the classes of gauge fixings under which it holds, clarifying the role of ghost fields in the 
formalism. 
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